Astronomy & Astrophysics manuscript no. ms 


© ESO 2010 


September 7, 2010 





On the stability of the thermal Comptonization index in 
neutron star low-mass X-ray binaries in their different spectral 

states 



o : 



R. Farinelli 1 & L. Titarchuk 1 ' 2 



Dipartimento di Fisica, Universita di Ferrara, Via Saragat 1, 44100 Ferrara, Italy 
p , 2 NAS A/GSFC, Greenbelt, MD 2077 1 

m : Preprint online version: September 7, 2010 



X 

Or 
i ' 

o ■ 

5h ■ 
-i— > . 

a . 



ABSTRACT 

Context. Most of the spectra of neutron star low mass X-ray binaries (NS LMXBs), being them persistent or transient, are characterized by 
the presence of a strong thermal Comptonization bump, thought to originate in the transition layer (TL) between the accretion disk and the NS 
surface. The observable quantities which characterize this component dominating the emission below 30 keV, are the spectral index a and the 
rollover energy, both related to the electron temperature and optical depth of the plasma. 

Aims. Starting from observational results on a sample of NS LMXBs in different spectral states, we formulate the problem of X-ray spectral 
formation in the TL of these sources. We predict a stability of the thermal Comptonization spectral index in different spectral states if the 
energy release in the TL is much higher than the intercepted flux coming from the accretion disk. 

Methods. We use an equation for the energy balance and the radiative transfer diffusion equation for a slab geometry in the TL, to derive 
a formula for the thermal Comptonization index a. We show that in this approximation the TL electron temperature kT c and optical depth 
to can be written as a function of the energy flux from the disk intercepted by the corona (TL) and that in the corona itself, Qdisk/Qcor- As 
spectral index a depends on kT a and to, this in turn leads to a relation a = /(Qdisk/Qcor), with a ~ 1 when Qdisk/Qcor >C 1. 
Results. We show that the observed spectral index a for the sample of sources here considered lies in a belt around 1 ± 0.2 a part for the 
case of GX 354—0. Comparing our theoretical predictions with observations, we claim that this result, which is consistent with the condition 
Qdisk/Qcor <C 1, can give us constraints on the accretion geometry of these systems, an issue that seems difficult to be solved using only the 
spectral analysis method. 
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It is well known since the 80's that the spectra of low-mass 
X-ray binaries (LMXBs) hosting a neutron star (NS) can be 
described up to about 30 keV by a two-component model rep- 
resenting the contribution from different emitting regions of 
the system. However, the interpretation of the spectra is not 
unique as demonstrated using the variety of different models 
used across the years. 

Before the BeppoSAX and RXTE era, two concurring mod- 
els were classically used to described X-ray emission in 
LMXBs. In the so-ca lled "eastern model" dMitsuda et alJl984r. 
Mitsud a et al. 1989) the spectra were fitted by the sum of a soft 



blackbody (BB) emission (actually modeled by a multi-colour 
disk BB spectrum) attributed to the accretion disk, plus a hot- 
ter simple or Comptonized BB claimed to originate close to the 
NS surface. On the other hand, in the "western model" inter- 



pretation dWhite et alJll986l:IWhite et al. II 19881) . the direct BB 
component was attributed to the NS surface, while an unsatu- 
rated Comptonization spectrum was thought to originate from a 
hot corona above the inner accretion disk, which supplies most 
of the soft seed photons for Comptonization. 

In fact, even after the advent of BeppoSAX and RXTE, 
the persistent emission of NS LMXBs was described by 
the sum of a BB component plus a thermal Comptonization 
(TC) sp ectrum, usually described by the XSPEC Comptt 



model (ITitarchukl 11994b iHua & Titarchuk I [1995). Despite 



the significant improvement in our knowledge of the source 
spectral properties by means of the broad-band observa- 
tions, the BB+TC model resulted and subjected to a di- 
chotomy. Indeed, both of the cases where the temperature of 
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the dir ec t BB s pectrum fcTbb is lower ( e.g.. iDi Salvo et al. 
2000allbL EoOlt foosterbroek et al. I l2Q0lL iGierlinski & Done 
200 2i;lLavagetto et al. ll2004IPaizis et al 112005b or higher (e.g. , 



Paizis et al. 1 120051: iFarinelli et al. I l2007r, iFarineili et al J 12008 
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hereafter F08) than that of the thermally Comptonized seed 
photons, kT s , provide in general equally acceptable good fits. 

The consequences of these results were the interpre- 
tation of the BB emission as due either to the accre- 
tion disk (fcTbb < kT s ) or to the NS surface (kT hh > 
kT s ). Moreover, in addition to the pers i stent X-ray emission , 
Beppo SAX (e.g.. bi Salvo et al. lEoOObt bi Salvo et al.1 l2002h 
RXTE dD'Amico et al. I l200ll bi Salvo et al.1 120061) and later 
also INTEGRAL dPaizis et al. 120061. hereafter P06) allowed the 
possibility to discover in bright NS LMXBs a transient power- 
law (PL) X-ray emission above 30 keV . 

Motivated by the need to put some order and give an uni- 
fied scenario in the different NS systems spectral states, Paizis 
et al. (2006, hereafter P06) performed a systematic observa- 
tional campaign with the ISGRI (20-200 keV) monitor onbo- 
rad INTEGRAL. Using long-term average s pectra and includin g 
former INTEGRAL results on GX 354-0 dFalanga et al. 112006b . 
P06 classified NS LMXBs into four main states: the hard/PL, 
low/hard, intermediate and soft. The high-energy (> 20 keV) 
spectra in different sources were interpreted by P06 as the re- 
sult of the interplay between thermal and bulk Comptonization 
processes whose relative efficiency is ultimately dictated by 
the mass accretion rate. At high energies (where the direct BB 
component is negligible), the hard/PL state spectra can be fit- 
ted with a simple PL component; the low/hard spectra by a TC 
spectrum of soft (^ 1 keV) BB-like photons off an electron 
population with kT c ^20 — 30 keV and tq 3; the interme- 
diate state spectra by TC spectrum with kT c ^3 — 5 keV and 
To ^ 5 plus a PL component with photon index r ~ 2 — 3; the 
soft state spectra by a TC component similar to the intermedi- 
ate state but without the high-energy X-ray tail. 

From the observational point of view, the quantities di- 
rectly measurable in the data are the cut-off energy E c of 
the dominating TC bump and the spectral slope, parametrized 
through the energy index a (=T — 1). For pure TC spectra, a is 
tightly co rrelated with the plasma temperature kT e and optical 
depth to (ISunvaev & Titarchuk 1980; iTitarchuk & Lvubarskiil 



hereafter TL95), while in the presence of a converging 
flow (bulk motion), the shape of the velocity field also de 



termines the emerging spe ctral slope (Titarchuk et al. 1 11997 



Laurent & Titarchuk 



19991 F08). In mathematical terms 



represents the index of the system Green's function (GF), 
namely as it responds to monochromatic line injection. The 
resulting emerging spectrum is then obtained as a convolu- 
tion of the GF with the input seed photon spectrum. We also 
emphasize that the hard/PL, hard and soft states spectra of 
NS LMXBs are characterized by the presence of just one 
Comptonization index (see Fig. 4 in P06). In the first case, its 
value is interpreted as a result of a mixed TC plus BC pro- 
cess, while in the latter two cases, a is derived from pure TC. 
On the other hand, the intermediate state spectra show two 
Comptonization indexes, one is related to the persistent TC 
component and the other one characterizes the transient PL- 
like hard X-ray emission. 

In Section|2]we give an overview of the theoretical and ob- 
servational issues related to spectral evolution in X-ray binary 
systems hosting a NS or a black hole (BH) outlining the dif- 
ferences among the two classes of sources. We subsequently 



report on results related to the observed TC index a for a sam- 
ple of NS sources. In Section[3]we propose a theoretical model 
based on diffusion formalism for radiative transfer in order to 
explain the observational results. In Section [4] we discuss the 
comparison between theory and data, while in Section [5] we 
draw our conclusions and give future observational perspects. 

2. Spectral index evolution in NS systems 

Starting from the considerations of the previous section, it is 
important to make a comparison between the index evolution in 
accreting BH and NS sources. The spectral state of BH sources 
may be generally divided into low/hard state (LHS), where the 
spectrum is dominated by a TC component with electron tem- 
perature kT c ~ 60-100 keV, intermediate state (IS), with a BB 
bump (presumably coming from the accretion disk) and a su- 
perposed PL high-energy component, and high/soft state (HSS) 
where the BB component gets even stronger and the PL emis- 
sion gets steeper. 

Actually, in BH systems generally one high energy photon 
index V is observed [see e.g. rece nt results on GRS 1915+015 



by Tit archuk & Seifina (2009) and Shaposhnikov & Titarchuk 
(2009)] and its value evolves with the source mass accre- 
tion rate. The latter, which is not a direct observable quan- 
tity, may be inferred in an indirect way by means of the the 
normalization of the disk flux A^isk- In a T vs Ad; s k/Quasi 
Periodic Oscillation (QPO) frequency diagram, it was found 
(SharJoshnikov& Titarchuk 2009; Titarchuk & Seifina 2009 



iMontanari et al. 1120091) that the photon index T progressively 
increases from ~ 1.6-1.8 as the source moves from the 
low/hard to the high/soft states, until it reaches saturation 
around T ~ 2.2 — 2.4 (depending on the source). This inter- 
pretation of the observed BH sp ectra have been perfo rmed us- 
ing "so called" the BMC model ( Titarchuk et al. II 1 997h . whose 
emerging spectral shape is described by the formula 



F(E) = 



Cn 



[BB{E)+AxG(E,E )*BB(E )}. (1) 



A+V 

In formula (Q}, the first term of the right-hand side represents 
the direct seed photon BB-like spectrum, while the second 
one gives its modification due to Comptonization (convolution 
with GF). The hardness of the spectrum is dictated by the en- 
ergy index a of the GF which, in BMC, is a broken-PL with 
G(E, E Q ) oc E a+3 for E < E„ and G{E,E ) oc E~ a for 
E > Eq, where Eq is the monochromatic input energy. It is 
important to keep in mind that the GF in BMC does not con- 
tain the exponential spectral rollover, so when the latter is ob- 
served in the X-ray spectrum of a source, it can be taken into ac- 
count by multiplying BMC with an exponential e-folding factor 
oc e~ E / B °. 

The index saturation observed in BH sources as they move 
to the HSS finds a natural explanation in the framework of a 
bulk-dominat ed scenario (see mathema tical proof of this state- 
ment in e.g. 



Titarchuk & Seifina 2009). Moreover, it is inter- 



esting to note that when the sources are in the LHS, the points 
in the T vs Adi S k/QPO diagram form a plateau around T ~ 1 .5 
before the rising phase. This clear mapping of the energy (or 
photon) index evolution is possible in BH sources as they are 
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generally strongly variable and their X-ray spectrum allows a 
one-to-one correspondence between the spectral state and the 
energy index a (or T) using spectral modeling according to for- 
mula ([T). One of the advantages of using BMC model is that it 
is a generic Comptonization model, namely it allows to map the 
spectral evolution of sources through the Comptonization index 
a, which is a direct measurable quantity, no matter which are 
the underlying physical conditions. The theoretical interpreta- 
tion of the source spectral formation is in fact related to a later 
scientific discussion. 

For NS sources however, situation is less straightforward 
to face. In fact, most variable sources such as, e.g., GX 354-0 
dFalanga et al. I 120061) or 4U 1608-52 jGierlinski & Done 



2002)) when moving from the hard state to the soft state do 



exhibit pure TC spectra, with electron temperature kT c pro- 
gressively decreasing and optical depth tq increasing, respec- 
tively. On the other hand, bright LMXBs of the GX class such 
as the classical six known Z sources (Sco X— 1, GX 17+2, Cyg 
X-2, GX 340+0, GX 5-1 and GX 349+2) and, more recently, 
GX 13+1, show only a small evolution of their persistent X-ray 
continuum (dominated by the strong TC bump with kT c ~ 3- 
5 keV and tq 3> 1). In addition they are characterized by the 
presence of a transient hard X-ray PL-like component (inter- 
mediate state). 

Other persistently bright sources such as GX 3+1, GX 9+1 
and GX 9+1 have been only observed in the soft state with no 
evidence of hard X-ray tails. In fact, mapping the evolution of 
the transient hard X-ray tail of NS sources in the intermediate 
state has not been yet possible because of the in sufficient statis 



importance. The limited broad-band resolution allows how- 
ever the possibility to have different models with same good 
fitting results. For example, recent results of the analysis of 
PC AJRXTE data for 4U 1608-52 by Ding et al. (2010 in prepa- 
ration) show that the soft state spectrum of the source can be 
fitted either by a two-BB model or by a single TC model (e.g., 
COMPTB) with a ~ 1. Note also that Falanga et al. (2006) fit- 
ted the soft state spectrum of GX 354-0 with a DBB+COMPTT 
model, from which the inferred a-value is again ~ 1 . Actually, 
when looking at the BeppoSAX results obtained over years of 
observations of NS LMXB sources, it results that the soft state 
spectra of these systems, rather than 2 BBs, need to be de- 
scribed by the sum of a BB component plus an unsaturated TC 
spectrum with cut-off energy below 10 keV. 

2.1. Observational results 

We considered a sample of sources taken from the literature 
and for which we can make a fiducial measurement of a. In 
our choice of the sources, we adopted the criterion to consider 
those in which a was determined either directly from the fit, 
as it can be done using the COMPTB model (where a is a free 
parameter, see F08), or can be derived from the te mperature 



and o ptical depth obtained by the COMPTT model (Titarchuk 

Il994h . 

The sources belonging to the first case are ScoX-1, 
GX17+2, Cyg X-2, GX 340+0, GX3+1 (se e Table 2 in 
F08) and GS 1826- 238 (see Table 2 i n ICocchi et al 



tics available at high energies. So me attempts C Di Salvo et al. 



bOOOatlDi Salvo et al.ll2002ll2006b were actually undertaken to 
establish changes in the intensity of the transient hard tail by 
fitting it with a simple PL and, when statistics was poor, fix- 
ing the index allowing to vary only its flux. But nothing could 
be concluded about the index evolution. Thus, at the present 
status of knowledge, serious investigations can be performed 
only on the evolution of the spectral index related to the per- 
sistent TC component. This was done for the first time by 



1201(1 For GX 349+2 dDi Salvo et al. lEoOOah and GX 354-0 
(Ei Salvo et al. 2000b) we derived the value of spectral index 
a using equation for the non-relativistic regime (see Eq. [22] in 
Titarchuk & Lyubarskij 1995, hereafter TL95) 



L 

e' 



(2) 



where = kT e /m e c 2 and /3-parameter defined in formula 
(17) of TL95 for spherical geometr y as it was assumed by the 



Titarchuk & Shaposhnikov I (|2005, hereafter TS05), who per- authors. In the case of X1658-298.IOosterbroek et al. I (120011) 



formed a systematic analysis of the variable NS X-ray binary 
4U 1728-34 from the hard state to the soft state using obser- 
vations from the Proportional Counter Array (PCA, 3-30 keV) 
onboard RXTE. The source spectra were fitted with a two-BMC 
model, with A ^> 1 (see formula (T|) and with the GF spectral 
index fixed equal for both BMC components. Thus, the model 
used by TS05 actually was 

F{E) = C ns G(E,E )*BB(E )+C disk G(E,E )*BB(E ). 

The main result found by TS05 was that as the source 
moved from the hard to the soft state the index T (= a + 1) pro- 
gressively increased with no evidence of saturation, unlike the 
BH case. Actually, the final soft state of 4U 1728-34was repre- 
sented by the sum of two BB component, because for a>l, 
G(E, E Q ) * BB(E ) » BB(E). 

This result however needs a revision. Namely, fitting the NS 
LMXBs soft state spectrum with a two-BB model can be pos- 
sible because of the lack of data below 3 keV, which is of key 



assumed a slab geometry, thus a was obtained from equation 
(f2]l an d formula ( 1 7 ) of TL 95 for a slab geometry. For IE 1724- 
3045, iBarret et al~l j2000l) report the best-fit value of optical 
depth To of Comptonization region both for the case of spher- 
ical and slab geometry. We checked that the two derived val- 
ues of a are perfectly consistent. The errors on a for sources 
for which kT c and tq were reported, have been computed con- 
sidering that the function a[kT e , f3(ro)] gets its absolute mini- 
mum and maximum values at the boundary of the box of its do- 
main delimited by the minimum and maximum value (fcT e mm , 
fcT" lax ) and (r™ ln , r™ ax ) obtained in computing the errors at 
90% confidence level for the electron temperature and optical 
depth by XSPEC. 

In Figure[T]we report the measured values of a for this sam- 
ple of sources as a function of the electron temperature kT e . 
This parameter can be considered indeed a good tracer of the 
source spectral state because kT c decreases when sources move 
from the hard to soft state as a result of a more efficient elec- 
tron cooling by the enhanced seed photon supply. Moreover, 
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Fig. 1. Thermal Comptonization index a for sources in dif- 
ferent spectral states as a function of the electron tempera- 
ture kT c . Reference papers: Cyg X-2, Farinelli et al. (2009) 
Sco X-1, GX 17+2, GX 340+0 and GX3+1, Farinelli et al 
(2008); GX 354-0, Di Salvo et al. (2000b); GX 349+2, Di 
Salvo et al. (2001), X 1658-298, Oosterbroek et al. (2001) 
GS 1826-238, Cocchi et al. (2010); IE 1724-3045, Barret et 
al. (2000). 



the electron temperature kT c is a directly measurable quantity 
because it is related to the cut-off energy of the spectrum and it 
presents the advantage of being distance-independent. On the 
other hand, the instrumental energy-band coverage and accu- 
mulation time can p lay some role in biasin g the measured index 
value. For instance. iFalanga et al. 1 d2006h performed a system- 
atic analysis of GX 354-0 as a function of its position on the 
hardness-intensity diagram when the source moved from the 
hard to soft state. The authors fitte d the 3-100 keV spe ctrum 



with a multicolour disk BB (DBB, iMitsuda et"ai]|l984l) pi 



COMPTT model using the slab geometry. We computed the de- 
rived a-values from their best-fit parameters but the trend was 
not monotonic, covering the range ^1-3 and reflecting such 
a behaviour of the kT c — To parameters (decreasing of kT e 
was not followed by increasing of tq, as expected). It is not 
clear whether this was due to the lack of data below 3 keV, 
which is necessary to constrain the seed photon temperature or 
due to accumulation time. Thus we prefer to skip these mea- 



surem ents. It is also worth mentioning that Oosterbro ek et al 



(2001) performed a BeppoSAX analysis of GX 3+1 and Ser 
X-1, both sources characterized by typical soft state spectra, 
fitting them with a DBB+COMPTT model, but they did not 
specify which geometry (sphere or slab) it was assumed. We 
found a sph = 2.44±° ^ and a sl = 0.87±° H for GX 3+1, 

0.20 



a sph = 1.411^;^ and a ls = 0.451^ for Ser X-1, respec 



tively. 

Looking at Fig. Q] we note that for all the analyzed sources 
the spectral index a lies in a belt around 1 ± 0.2, a part for the 
case of GX 354-0 where a ~ 1.6. We give a possible interpre- 
tation of these observational results in the Discussion. 



3. A model of Comptonization region in the 
neutron star case 

The determination of spectral index a obtained from 
Comptonization of seed photons in a bounded medium has 
been faced since a long time. The emerging radiation spec- 
trum depends on several parameters such as the geometry of 
the plasma (e.g., slab or sphere), the electron temperature and 
optical depth, and the space distribution of the seed photons in- 
side the medium. Sunyaev & Titarchuk (1980) report the value 
a obtained from the solution of the stationary radiative trans- 
fer equation in the non-relativistic case (Fokker-Planxk ap- 
proximation) obtained as a convolution of the time-dependent 
equation with the time-escape probability distribution P(u) for 
the case of source photons distributed ac cording to t he firs t 
eigenfunction of t he spa ce operator. Later iTitarchuk I (1 19941) . 
Hua & Titarchuk I (11995b and TL95 extended the results to the 
sub-relativistic case considering both slab and spherical geom- 
etry. 

In order to understand what it does happen in NS LMXBs 
sources, one has to consider the hydrodynamical conditions in 
the region between the accretion disk and the NS surface. We 
will refer to this region as the transition layer (TL), often also 
called boundary layer. Actually, the production of a strong TC 
bump in the persistent X-ray spectra of NS LMXBs is thought 
to originate in this TL, namely the region where matter deviates 
from its Keplerian angular velocity in order to match that of the 
slowly spinning NS. The radiative and hydrodynamical config- 
uration of the TL is mostly dictated by the Reynolds nu mber 7 
dTitarchuk et al. Ill 9981; ITitarchuk & Osherovich II 19991) . which 
is proportional to the mass accretion rate and is eventually the 
inverse of the viscosity parameter of the Shakura-Sunyaev disk. 
In particular, 7 determines the radial extension of the TL and in 
turn, from the mass accretion rate, its optical depth. It is worth 
pointing out that determination of the vertical height-scale of 
the TL is in fact a very complicated problem, as it should re- 
quire a complete 3D magneto-hydrodynamical treatment of the 
problem. Using the slim disk (thus vertical -averaged) equa- 
tions for determining the radial thermo-h ydrodynamical struc 



tureo f the TL may be in fact an issue (e.g. JPopham & Sunyaev 
20011). 



The enhanced radiation and thermal pressure because of 
higher electron temperature are expected to increase the ver- 
tical height-scale of the TL with H « R ns . Moreover, the 



solution for the angular momen t um eq uation (ITitarchuk et al 



1998; ITitarchuk & Osherovich I H 999) for Reynolds number 
7 > 5 — 10 gives a TL radial extension ARtl Si 0-5Rns ■ 
With these characteristic length-scales, it seems more plausible 
to approximate the TL geometry to a slab whose normal is di- 
rected along the disk plane. It is worth emphasizing that Haardt 
& Maraschi (1993, hereafter HM93) determined the theoret- 
ical Comptonization index a considering a two-phase model 
for accretion disks in AGN in which a hot corona is surround- 
ing and sandwiching the underlying cold accretion disk. The 
model can in principle be applied also to the case of solar-mass 
BH sources. HM93 assume that the corona and the disk are two 
slabs at significant different temperatures and put in contact 
each other. The authors concentrated on the case of high tem- 



Farinelli & Titarchuk: Thermal Comptonization index in NS LMXBs 



5 



perature (kT c <; 50 keV) and low optical depth (t <1) for the 
corona, so that the diffusion approximation cannot hold, unlike 
what we are considering. One of the consequences of the high- 
temperature treatment is that electron scattering is anisotropic 
with a significant fraction of the power back-irradiating the 
disk. In the HM93 model, the inner boundary condition of 
the hot corona is the disk cool surface (with fclbb < 5 eV) 
with energy-dependent albedo. Note also that in their geome- 
try, 100% of the disk flux is intercepted and reprocessed by the 
top plasma. In the geometry considered here on the other hand, 
it is possible that part of the disk emission directly escapes the 
system, while a fraction of its flux is intercepted by the TL. We 
are actually interested here on this portion of the intercepted 
disk flux (see next Section). 

Because of these differences between the two models, a di- 
rect comparison of the derived theoretical results is not straight- 
forward. The reader can refer to the paper of HM93 for further 
details in order to better understand the differences of our and 
their approaches. 

3. 1. Energy release in the NS transition layer 

The energy balance in the TL is dictated by Coulomb colli- 
sions with protons (gravitational energy release), while inverse 
Compton and free-free emission are the main cooling chan- 
nels (see a formulation of this problem in the pioneer work 
by Zel'dovich & Shakura 1969). In fact, for the characteristics 
electron temperature (3 keV Ss kT e 30 keV) and density val- 
ues (Ss 10~ 5 g cm' 3 ) of these regions for LMXBs, Compton 
cooling dominates over free-free emission and the relation be- 
tween the energy flux per unit surface area of corona Q COI , the 
radiation energy density e(r) and electron temperature T e is 
given by (see also TLM98) 



Qc 



TO 



20.2e(r)T e (r) 



(3) 



where tq is the characteristic optical depth of the TL. 

The distribution e(r) is obtained as a solution of the diffu 
sion equation 

dre = ggtgt 
dr 2 ctq ' 

where now Q to t = Qcor + Qdisk is the sum of the corona 
(TL) and intercepted disk fluxes, respectively. The two bound- 
ary conditions for equation (0]) are written as 



(4) 



de 



0. 



de 3 

Tt 2 £|t=0 







(5) 



(6) 



which represent the case of albedo A=l at the NS source 
(r = To) and no diffusion emission falling from outside onto 
outer corona boundary (t = 0). The condition for A=l arises 
from the well-established observational result of NS tempera- 
ture fcTbb ~ 1 keV, which implies the presence of a ionized 
NS atmosphere. This is different from the case considered by 
HM93, where the cool disk temperature (< 5 eV) gives rise to 
an energy-dependent albedo with photoelectric absorption for 



impinging photons with energy 5s 10 keV. Another important 
consideration to keep in mind is that the diffusion equation dU 
is to be considered frequency-integrated. This means that we 
are not dealing with the specific (energy-dependent) shape of 
the reflected spectrum from the NS surface, we are only deal- 
ing with the total energy density. 

The solution for e(r) is then given by 



e(r) 



2Q t 



1 



-TO 



Tl) 



(7) 



It is worth pointing out that de/dr > for r < To, and as 
F ra d oc de/dr for NS sources the radiative force always plays 
against gravity, unlike the case of BH sources. 

Note that the spectra of NS sources both in the soft 
and hard state can be adequately fitted by single-temperature 
Comptonization models (Pazis et al. 2006, Falanga et al. 2005, 
Farinelli et al. 2008, Cocchi et al. 2010). This observational fact 
demonstrates that an assumption of isothermal plasma in the 
TL can be applicable to X-ray data analysis from NS binaries. 
The question is how one can estimate this average temperature 
of the TL which is, in fact, established by photon scattering and 
cooling processes. 

In order to establish this average plasma temperature T e one 
should estimate the mean energy density in the TL as 



T Jo 



:{r)dT 



-(2 + r ). 



(8) 



It is worth emphasizing the similarity between equa- 
tions 01) and (EH) in our p aper and equation (13) in 



Bisnovatvi-Kogan et al. 



1980) who studied the radiation 
emission due to gas accretion onto a NS. If we now substi- 
tute the result of equation ([8]) into equation (0), after a bit of 
straightforward algebra we obtain 



kT eTo {2 + To) 



0.25 



m e c* 



l + Qdisk/Q c 



(9) 



Keeping in mind the definition of the Compton param- 
eter Y w AN SC (Rybicki & Ligthman 1989), where A ~ 
AkT e /m e c 2 and N sc ~ Max(ro,To) are the average photon 
energy gain per scattering and average number of scatterings, 
respectively, we can rewrite equation (01 as follows 



Y 



1 



(10) 



1 + Qdisk/Q cor 

Equation ( TTOb is one of the main points of our theoreti- 
cal model and shows that in the diffusion approximation the 
Compton parameter, which determines the spectral index, is 
just a function of the corona and disk cooling fluxes. 

3.2. The radiative transfer model of small variations 
of index a during spectral state transition in NS 
sources 

As we have shown in Section 2.1, the observed spectral in- 
dex a of most NS LMXBs undergoes small variation around 1, 
namely a = 1 ±0.2 when the electron temperature of Compton 
cloud varies from about 2.5 to 25 keV (see Fig. Q]). Thus here 



6 



Farinelli & Titarchuk: Thermal Comptonization index in NS LMXBs 



we propose a model the spectral formation in the TL (corona) 
which can explain the stability of index a if energy release in 
the disk is much less than that in the TL. Namely we show that 
a w 1 + 0(Qdisk/Q coiv ■ 

As already pointed out in classical works (Sunyaev & 
Titarchuk 1980, 1985, hereafter ST85, Titarchuk 1994), spec- 
tral formation in plasma clouds of finite dimensions (bounded 
medium) is related to the distribution law of the number of scat- 
terings that seed photons experience before escaping. If u av 
denotes the average number of photon scatterings and the di- 
mensionless scattering number is u = N e arct, then the distri- 
bution law for u 3> u av is given by (see ST85) 



P(u)=A(u av ,T )e-( )u . 



(11) 



For a diffusion regime when To <; 1.5, it results B — 
where Ai is the first eigenvalue of the diffusion space opera- 
tor. As reported in ST85, the eigenvalue problem for photon 
diffusion in a slab with total optical depth 2to is derived from 
solution of the differential equation for the zero-moment inten- 
sity 



with absorption boundary conditions dJ/dr — (3/2) J = and 
dJ/dr + (3/2)J = 0, for r = and r = 2to, respectively. 
This leads to the trascendental equation for the eigenvalue A n , 
n = 1,2,3... 

2 

tan(A n T ) = — , (13) 

6A n 

which has the solution for To 3> 1 and n = 1 

TT 



Ai = 



2(t + 2/3) ' 



(14) 



It is important to emphasize that the same result for Ai is 
obtained by solving equation (TT2l for a slab with total optical 
depth To but with reflection condition dJ/dr = at t = To. 
This is not surprising as such a condition is actually met at 
the center of a symmetric slab with total optical depth 2to and 
< t < 2to. Thus, the same mathematical result is obtained 
for two different geometrical configurations. In the first case 
(symmetric slab with total optical depth 2to) it represents, e.g., 
an accretion disk (ST85 treatment), in our present case we are 
dealing with a boundary layer with total optical depth To, whose 
asymmetry is due to the presence of a reflector (NS surface) at 
one of the two boundaries. In both cases one obtains 



12(t + 2/3)2" 



(15) 



Generalizing to the case of arbitrary optical depth To, the 
diffusion operator Ldiff = (l/3)d 2 J/<iT 2 is replaced by the 
radiative transfer operator L T applied to J(t) (see ST85 and 
TL95) which for the disk geometry is 



L T J=\ 



2r 



J(r')Si(|r-r'|)dr'- J, 



(16) 
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Fig. 2. Values of the /3-parameter as a function of the optical 
depth To. The solid, dashed and dotted lines correspond to defi- 
nition of B given in equations ( TTBI ), ( fTTI i and (1181 1. respectively. 



where E\ (z) is the exponential integral of the firs t order. In this 
case, the derived value for (B is (Titarchuk 1 11 994 TL95) 



(12) P 



-X.Z5t q \ 



a , 10 

0.45e- 3 7To ln . (17) 

3t 



12(t + 2/3) 2 

Now having in mind equation (0, we introduce the param 



eter 



/?diff 



1 



(18) 



r (2 + T )' 

and in Figure|2]we show the values of B for the cases reported 
in formulas (15[ , ( TTTb and (fT8l as a function of optical depth 
To- It is possible to see that actually for To <; 1.5 all values of 
B are practically close each other, but they deviate for To ^ 1. 
For example their difference is about 30% for To=l. 

Using the definition of a (see Eq.[2j, where 8 is replaced by 
Bdifi (Eq.QjOl, and equation (0, we obtain the diffusion spectral 
index as 



«diff 



l + Qdisk/Q c 
025 



(19) 



Of a d ifF = 1 + 0.8 Qdisk/<3cor for Qdisk/Qcor < 1. 

Thus as it follows from Eq. ( fT9l ), in the diffusion regime the 
TC spectral index can be expressed in terms of Q d isk/Qcor (the 
intercepted disk over coronal fluxes), instead of TL electron 
temperature kT e and optical depth To (see Eqs. |2] and fTBI ). In 
Figure|3]we present a plot of a d ifF as a function of Qdisk/Qcor, 
which shows that it ranges from 1 to 1.6 as Q d i S k/Qcor in- 
creases from to 1 . One can see the observable values of index 
a ~ 1 takes place if the energy release in the disk is much less 
that in TL, namely if Qdisk/Qcor < L 

4. Results and discussion 

In this Paper we compare observational data of a sample of 
NS sources (Fig. [T) with the theoretical results which follows 
from the radiative transfer model in the diffusion approxima- 
tion. The data show that independently of the source spectral 
state, which we have parametrized through the measured TL 
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Fig. 3. Theoretical thermal Comptonization index a as a func- 
tion of the ratio Qdisk/Qcor according to equation ( fT9l ). 



electron temperature kT e , the spectral index a = 1 ± 0.2. We 
derived an estimate of the energy index a for TC spectra in NS 
LMXBs using an equation for the diffusion approximation in a 
slab geometry with the reflection (100 %) boundary condition, 
which is valid for optical depth tq > 1.5. In particular, we find 
that in this approximation it is possible to express the value of 
a as a function of the ratio of the flux from the accretion disk 
intercepted by the corona (TL, see Eq. fl9l ) and the energy 
release in the corona itself. The agreement of the model with 
the data takes place when the condition Qdisk/Qcor "C 1 holds 
(see Fig.0. 

This behavior of spectral index a actually has important 
consequences on putting constraints in the accretion geom- 
etry of NS LMXBs. In fact, as already pointed-out in the 
Introduction, the broad-band persistent X-ray spectra of these 
sources are usually fitted by a two-component model consist- 
ing of BB-like emission plus a strong TC bump. Both the cases 
for which the BB temperature (fcTbb) is lower or higher than 
that of the seed photons of the TC bump (kT s ) provide equally 
acceptable fits. In the first case (fcTbb < kT s ), the origin of the 
direct BB component is attributed to the accretion disk, in the 
second case (fclbb > kT s ) to the NS surface. This dichotomy 
in the spectral analysis has not been overcome for a long time. 
Some help in this direction has come with the discovery of the 
transient X-ray tails in some of the brightest sources (see refer- 
ences in the Introduction): in fact, if the origin of the hard PL- 
like X-ray emission is attributed to a combined thermal plus 
bulk (converging flow) effect in the innermost part of the TL 
close to the NS surface, then it becomes natural to suggest that 
the observed direct BB photon spectrum mostly originates in 
the TL/NS surface region, pr oviding the seed photo ns for bulk 
Comptonization (see Fig. 1 in Farinelli et al. I l2007h . 

The theoretical results derived in this Paper, in our opinion 
strengthen this scenario. In fact, when computing the total en- 
ergetic spectral budget of the sources in the 0.1-40 keV (where 
most of the emission is produced) the dominating TC bump car- 
ries out more than 70% of the source luminosity, the remaining 
part due to the direct BB component. If this BB originates close 
to the NS surface, then it is evident that the disk contribution to 



the X-ray luminosity is very small. Given that Qdisk in equation 
([T9l > represents the flux from the accretion disk intercepted by 
the corona and thus is smaller or at most equal to the directly 
emitted part, this eventually leads toa~ 1. In this framework, 
we can make some considerations about the higher value of a 
(~ 1.6) measured for GX 354-0 (see Fig. Q}. Di Salvo et al. 
(2000b) fitted the broad-band spectrum of the source with a 
BB plus COMPTT model (Titarchuk 1994) with kT hh < kT s , a 
modelization corresponding to the case in which a significantly 
higher fraction of the X-ray luminosity comes from the accre- 
tion disk. This would turn of course in an enhanced value of 
Qdisk, and looking at Figure[3]it is straightforward to see that 
increasing Qdisk/Qcor corresponds to increasing a. Actually, 
it would be interesting to see what it does happen by fitting 
the BeppoSAX spectra of GX 354-0 with the same model but 
with kT hh > kT s . Note also in Fig. Q] that for IE 1724-3045 
and two spectral states of GS 1826-238, a ~1 with electron 
temperature kT c i; 20 keV. Using equation (|9]i and the best-fit 
values of kT c reported for the two sources, with Qdisk/Qcor ~ 
0, we obtain r ~ 1.3 for IE 1724-3045 and t ~ 1.6-1.7 for 
GS 1826-238, respectively. These values of the optical depth 
allow actually to deal with the diffusion approximation within 
degree of accuracy still satisfactory as can be seen from Fig. 
|2] Moreover, it is worth emphasizing that for a slab with opti- 
cal depth To and inner reflection boundary condition, photons 
which are back-scattered from the reflecting surface before es- 
caping experience in fact an optical depth ~ 2tq, which further 
enhances the diffusion approximation validity. 

The other issue to point out is that the lower limit on a (for 
the extreme case Qdisk/Qcor =0) derived by our model is 1, 
while there is a handful of sources for which 0.8 < a < 1. 
Different reasons can in fact concur to this result. First of all, it 
is well known that multi-component modeling of X-ray spec- 
tra may have some influence in the determination of the best-fit 
parameters. In particular, in the energy band where TC domi- 
nates (^ 30 keV) Gaussian emission lines around 6.4-6.7 keV 
are often observed and inclusion of narrow-feature component 
in the model can affect the continuum parameters. Additional 
biases can come from the energy-band coverage (in particular 
when using RXTE data, which start from about 3 keV) and un- 
certainties in the calibration of the instrumental effective area, 
which may play some role in particular for spectra which are 
far away from being powerlaw-like. 

In terms of theoretical predictions, our analytical model is 
intended to provide a description of the observed stability of 
the spectral index but may of course have some limitations, 
in particular the energy-independent treatment of the radiation 
field and the simple slab approximation for the TL geometry. 
However it is possible that the vertical height-scale H of the TL 
is higher than its radial extension, it is also likely that there is 
some dependence of H on the radial distance from the NS sur- 
face. Moreover, in the pure plane-parallel geometry photons are 
allowed to escape only through the surface of the slab, while in 
this case the slab has limited extension and photons presumably 
can escape also from its "walls". 
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5. Conclusions 

We reported results on the value of the thermal Comptonization 
(TC) spectral index a for a sample of LMXBs NS sources in 
different spectral cases and we found that a part from GX 354- 
it lies in a belt around 1 ± 0.2. We proposed a simple theoret- 
ical model using the diffusion approximation where a is found 
to be a function only of the ratio of the disk and corona fluxes. 
In particular when Qdisk/Qcor <C 1 it results a <~1, which 
is consistent with observations. We are hopeful that our work 
will encourage to significantly extend the sample of observed 
NS sources using archival data and observations from present 
and future missions, in particular using as broad as possible en- 
ergy band, especially below 3 keV in order to avoid biases in 
the spectral results. We also claim that our model can be helpful 
in solving the dichotomy related to the fact that equally good 
fit are obtained for the cases in which the observed direct BB 
component has a temperature higher or lower than that of the 
seed photons subjected to TC. 
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